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Abstract 

The q — 10 and q = 200 state Potts models coupled to 2d grav- 
ity are investigated numerically and shown to have continuous phase 
transitions, contrary to their behavior on a regular lattice. Critical ex- 
ponents are extracted and possible critical behavior for general g-state 
Potts models coupled to gravity is discussed. 



1 Introduction 

In the continuum formulation of 2d quantum gravity one is interested in 
studying the coupling of conformally invariant matter to gravity. The con- 
formal matter is characterized by its central charge c but yet it has only 
been possible to solve theories with c < 1 analytically. 

This has triggered a lot of numerical work on discretized models of 2d 
gravity with c > 1 [|l|. One such discretization method is dynamical trian- 
gulations where the integration over space-time manifolds is replaced by a 
sum over piecewise linear surfaces constructed by gluing together equilat- 
eral triangles. The matter is simulated by spins, living either on triangles 
or vertices. The spin models usually used are the q = 2, 3 and 4-state Potts 
models and the Gaussian model. All these models have continuous phase 
transitions on regular lattices, with divergent correlation length, and hence 
can be associated with some conformal field theories in the continuum limit. 
The above mentioned models all correspond to conformal theories with c < 1 
but taking multiple copies of them allows one to study numerically theories 
with c > 1 coupled to gravity. For unlike in the continuum formulation the 
discretized models are well defined for c > 1. 
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For g > 4 the Potts models have a discontinuous or 1st order phase 
transition on regular lattices and hence no corresponding conformal field 
theories. So until now coupling them to gravity has been discarded as unin- 
teresting. But a recent mean field solution of the q = oo state Potts model 
coupled to gravity showed the model to have a 3rd order magnetization 
transition from a high-temperature branched polymer phase to a low tem- 
perature pure gravity phase (c = 0) Moreover the model was shown to 
be equivalent to infinite copies of Ising models coupled to gravity (a c = oo 
model). 

By expanding around the q = oo solution this analysis was extended 
down to finite values of q. There two phase transitions were found with high- 
and low-temperature pure gravity phases separated by a branched polymer 
phase. The critical behavior of those transitions was identical, both had 
the specific heat exponent a = —1. As the value of q is decreased the two 
phase transitions merge at some critical value qc and the branched polymer 
phase disappears. At that value of q (estimated to be qc ~ 120 in Q) the 
perturbation expansion breaks down. 

It is worth comparing this with the situation for multiple q < 4 state 
Potts models coupled to gravity. For c < 1 the interaction between gravity 
and matter is weak, for all values of the coupling except f5c we have a pure 
gravity phase. But for larger values of c numerical simulations indicate 
a high- and low-temperature pure gravity phase with a branched polymer 
phase in between, i.e. the same scenario as for high-g. 

So the similarity between these two cases might indicate that the study 
of Potts models for large value of q might give us informations about mod- 
els with large central charge coupled to gravity. It might be possible that 
after coupling to gravity these model would describe some conformal matter 
coupled to gravity, the only problem being that we have not yet identified 
these conformal theories. 

This has prompted us to look at the g = 10 and q = 200 state Potts 
models coupled to gravity applying Monte Carlo simulations Using 
micro-canonical simulations (fixed area) the models are defined by 



where ai G {1, are the Potts spins, i is a lattice triangle, {ui} a spin 

configuration on the triangulation T and T(N) an appropriate class of tri- 
angulations (of size N). In these simulations we include triangulations al- 
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Figure 1: (a) The specific heat Cy versus /3, for different lattice sizes, for 
the q = W state Potts model coupled to gravity, (b) The energy cumulant 
Vl versus /? for the same model 

lowing sites joined by more than one link and sites connected to them selves 
as these can be shown to reduce finite size effects [^]. The triangulations 
were updated using the link-flip algorithm and the spin models with the 
Swendsen-Wang cluster algorithm. The lattice sizes used were between 250 
and 8000 triangles and, as the auto-correlations increase linearly with q, we 
used 5 X 10^ and 2 x 10^ sweeps for each /3 value for the q = W and q = 200 
state Potts model respectively. 

2 The order of the phase transitions 

The first interesting question is the order of the phase transitions, are they 
continuous or discontinuous after coupling to gravity. This we determine 
by looking at the scaling of the maximum of the specific heat with lattice 
size. If a phase transition is 1st order Cv{max) should scale linearly with 
the volume. Prom fig. la, where we plot the specific heat versus (3 for the 
q = \Q state Potts model, we see that the maximum increases slowly with 
lattice size. Looking at the quantity Cv{max))/N it approaches zero in the 
infinity volume limit, which is a strong indication for a continuous phase 
transition. 

Another quantity sensitive to the order of the phase transition is the 
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q PcjCy) PcjBC) 
10 1.412(1) 1.409(2) 
200 2.633(1) 2.629(3) 



"To 2.01(5) 2.03(4) 0.53(1) 1.11(2) 0.02(4) 0.02(3) 
200 2.49(15) 2.52(12) 1.18(4) 0.31(5) -0.83(6) -0.80(5) 



Table 1: Measured critical exponents for the q = 10 and q = 200 state Potts 
models coupled to gravity. For vdn we have values from (a) the scaling of 
m.Six{dBC/d(3} and (b) the scaling of M and dM/d(5. In a same way for a 
(a) comes from Cv{'rnax) and (h) from Cy(/3c)- 



energy cumulant 

If the transition is 1st order Vl should reach a non-zero minimum in the 

critical point, otherwise it is zero for all values of (3. In fig. 16 wc show that 
for the q = 10 state Potts model no indication of a non-zero minimum can 
be seen. 

Exactly the same is seen for the q = 200 state Potts model, and from 
these measurements we conclude that both models have continuous mag- 
netization transitions after coupling to gravity. It it might at first seem 
strange that the order of the phase transition can change that drastically 
with coupling to gravity, but we should remember that even for one Ising 
model (which can be solved explicitly) there is a change from a 2nd to a 3rd 
order transition with coupling to gravity. 



3 Critical exponents 

To locate the critical temperatures for the models we used the finite size 
scaling behavior of two quantities: (a) location of the peak in the specific 
heat and (h) the intersection of Binder's cumulant for different lattice sizes. 
Both quantities are expected to approach /3c as N'^^^^f^ . The fits to the 
scaling behavior are made easier as vdn can be determined directly from 
Binder's cumulant, as the maximum of its slope scales as N^'^" . The results 



4 



are shown in table 1. 

We then apphed standard finite size scahng in /?c to extract the the crit- 
ical exponents (3 (from the magnetization M), udn (from dM/df3), 7 (from 
the magnetic susceptibility x) a (both from the scaling of Cv{max) 
and Cy(/3c))- The results are shown in table 1. For g = 10 we tested that 
the scaling of the specific heat fits equally well to logarithmic divergence 
(Cy ~log(iV)). 

4 Discussion 

It is an interesting observation that the critical exponents for the q = 10 
state Potts model are very close to that of the g = 4 state Potts model 
coupled to gravity (which has (3 = 1/2, 7 = 1, a = and udn = 2). As 
the value g = 10 is chosen arbitrary this might indicate that there exists an 
whole range of q values, starting at g = 4 and ending at some qc, where we 
have the same critical behavior after coupling to gravity. 

On the other hand the g = 200 state Potts model has a close to — 1, 
so we might be in the high-g region were the calculation in Q is valid. 
Then we would expect a branched polymer phase and two phase transitions. 
Unfortunately we have only been able to locate one, but as the other one is 
not a magnetization transition it is not clear how to identify it. But for an 
interval of /?, before /3c, we see some evidence of a branched polymer phase, 
supporting the picture in Q . 

From this the following critical behavior of the g-state Potts models 
coupled to gravity sounds reasonable: For g = 4 up to some g^ we have the 
same critical behavior, i.e. a 2nd order phase transition (a = 0) where the 
details of the Potts models are less important than the coupling to gravity. 
And for g > gcj we also have the same critical behavior, this time with 
(presumable two) 3rd order phase transitions (a = —1) and a branched 
polymer phase. Whether gc^ = gcj or if we have some different behavior at 
intermediary values of g is not determined in these simulations. But it is 
clear that we need some theoretical understanding of what kind of conformal 
matter (if any) the g > 4 state Potts model coupled to gravity describes. 
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